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Abstract 

We give several equivalent combinatorial descriptions of the space of states for 
the box-ball systems, and connect certain partition functions for these models 
with the g-weight multiplicities of the tensor product of the fundamental rep- 
resentations of the Lie algebra 0l(n). As an application, we give an elementary 
proof of the special case t = 1 of the Haglund-Haiman-Loehr formula. Also, 
we propose a new class of combinatorial statistics that naturally generalize 
the so-called energy statistics. 
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1 Introduction 



The purpose of the present paper is two-fold. First of all, we would like to give an 
introduction to the beautiful Combinatorics related with Box-Ball Systems, and sec- 
ondly, to relate the latter with the "Classical Combinatorics" revolving around trans- 
portation matrices, tabloids, the Lascoux-Schiitzenberger statistics charge, Macdon- 
ald polynomials, [3l],[35], Haglund-Haiman-Loehr's formula [8], and so on. As a 
result of our investigations, we will prove that two statistics naturally appearing in 
the context of Box-Ball systems, namely energy function and tau-function, have nice 
combinatorial properties. More precisely, the statistics energy E is an example of 
a generalized machonian statistics [22], Section 2, whereas the statistics r related 
with Kostka-Macdonald polynomials, see Section 5.2 of the present paper. 

Box-Ball Systems (BBS for short) were invented by Takahashi-Satsuma [lH HI] 
as a wide class of discrete integrable soliton systems. In the simplest case, BBS 
are described by simple combinatorial procedures using box and balls. Despite its 
simple outlook, it is known that the BBS have various remarkably deep properties; 

• Time evolution of the BBS coincides with isomorphism of the crystal bases 
[T0| [5] . Thus the BBS possesses quantum integrability. 

• BBS are ultradiscrete (or tropical) limit of the usual soliton systems [IHl [29] . 
Thus the BBS possesses classical integrability at the same time. 

• Inverse scattering formalism of the BBS coincides with the rigged configuration 
bijection originating in completeness problem of the Bethe states [281 138] . 

Let us say a few words about the main results of our paper. 

• In the case of statistics tau, our main result can be formulated as a computation 
of the corresponding partition function for the BBS in terms of the values of 
the Kostka-Macdonald polynomials at t = 1. 

• In the case of the statistics energy, our result can be formulated as an inter- 
pretation of the corresponding partition function for the BBS as the g-weight 
multiplicity in the tensor product of the fundamental representations of the 
Lie algebra Qi{n). We expect that the same statement is valid for the BBS 
corresponding to the tensor product of rectangular representations. 

We are reminded that a q-analogue of the multiplicity of a highest weight A 
in the tensor product (S)a=i K„a;^„ of the highest weight Sa oj^^, a = 1, . . . , L, 
irreducible representations Vg^ujra algebra fl[(n) is defined as 

L 

(7-Mult [Vx : (g) = J2 ^v,R ^^vM^ 

a=l r) 

where K^^^r stands for the parabolic Kostka number corresponding to the se- 
quence of rectangles R := {(s^'')}a=i,...,L, see e.g. [22], [25] . 
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• We give several equivalent descriptions of paths which appear in the description 
of partition functions for BBS: in terms of transportation matrices, tabloids, 
plane partitions. We expect that such interpretations may be helpful for better 
understanding connections of the BBS and other integrable models such as 
melting crystals [3l], g-difference Toda lattices [6], ... . 

Our result about connections of the energy partition functions for BBS and q- 
weight multiplicities suggests a deep hidden connections between partition functions 
for the BBS and characters of the Demazure modules, solutions to the g-difference 
Toda equations, cf.[6]. 

As an interesting open problem we want to give raise a question about an in- 
terpretation of the sums Kri,R K^^x{Q,t), where Kn^\{q,t) denotes the Kostka- 
Macdonald polynomials [31], as refined partition Junctions ioi the BBS corresponding 
to the tensor product of rectangular representations R = {(s^")}i<a<„. See Conjec- 
ture 15.191 In other words, one can ask: what is a meaning of the second statistics 
(see [H]) in the Kashiwara theory [18] of crystal bases (of type A) ? 

Organization of the present paper is as follows. In Section 2, we review neces- 
sary facts from the Kirillov-Reshetikhin crystals. Especially we explain an explicit 
algorithm to compute the combinatorial R and the energy function. In section 3, we 
introduce several combinatorial descriptions of paths. Then we define several statis- 
tics on paths such as Haglund's statistics, energy statistics E and tau statistics t^''^. 
In Section 4, we collect necessary facts from the BBS which will be used in the next 
section. In Section 5, we present our main result (Theorem 15.71) as well as several 
relating conjectures. We conjecture that t^''^ gives independent statistics depending 
on one parameter r although they all give rise to the unique generating function up 
to constant shift of power. In Section 6, we show that the energy statistics E belong 
to the class of statistics r^''^ (Theorem 16.31) . Therefore r*"'* gives a natural extension 
of the energy statistics E. 

2 Kirillov-Reshetikhin crystal 

2.1 An^ type crystal 

Let Ws^"^ be a Uq{Q) Kirillov-Reshetikhin module, where we shall consider the case 
Q = An\ The module Ws'^^ is indexed by a Dynkin node r G / = {l,2,...,n} 
and s G Z>o- As a ?7g(A„) -module, Ws^^ is isomorphic to the irreducible module 
corresponding to the partition (s^). For arbitrary r and s, the module Ws^"^ is known 
to have crystal bases [HI [17], which we denote by S^'**. As the set, B^'^ is consisting 
of all column strict semi-standard Young tableaux of depth r and width s over the 
alphabet {1, 2, . . . , n + 1}. 

For the algebra An, let P be the weight lattice, {Aj G P\i G /} be the funda- 
mental roots, {ai G P\i G /} be the simple roots, and {hi G Homz(P, Z)|i G /} be 
the simple coroots. As a type An crystal, B = fi*"'* is equipped with the Kashiwara 



3 



operators Cj, : B — > B U {0} and wt : B — > P {i E I) satisfying 

f.{b) =b' ^ e,{b') =b if 6, 6' G B, 
wt im) = wt(fo) - a, if Mb) e B, 
{hi,wtib)) =^i{b)-ei{b). 

Here (■,■) is the natural pairing and we set ei{b) = max{m > | e"''b ^ 0} and 
= max{m > | ^ 0}. Actions of the Kashiwara operators Cj, for z G / 
coincide with the one described in [19]. Since we do not use exphcit forms of these 
operators, we omit the details. See [33] for complements of this section. Note that 
in our case An-, we have P = Z""*"^ and = — e^+i where is the i-th canonical 
unit vector of Z^^^. We also remark that wt(6) = (Ai, ■ ■ ■ , A^+i) is the weight of b, 
i.e., Aj counts the number of letters i contained in tableau b. 

For two crystals B and B', one can define the tensor product B ^ B' = {b ® b' \ 
b & B,b' ^ B'}. The actions of the Kashiwara operators on tensor product have 
simple form. Namely, the operators Cj, fi act on B ® B' by 

ei{b®b') 
fiib®b') 

and wt(& ® b') = wt(6) + wt(6'). We assume that ® 6' and b ® as 0. Then it is 
known that there is a unique crystal isomorphism R : B^'^ ® B'"''^' ^ B'^'''^' ® S'''*. 
We call this map (classical) combinatorial R and usually write the map R simply 
by ^. 

Let us consider the affinization of the crystal B. As the set, it is 

AS{B) = {b[d] \beB,deZ}. (1) 

There is also explicit algorithm for actions of the affine Kashiwara operators eo, 
/o in terms of the promotion operator [12]. For the tensor product b[d\ ® b'[d'] G 
Aff (S) (g) Aff (S'), we can lift the (classical) combinatorial R to affine follows: 

b[d] ® b'[d'] ~ b'[d' - H{b ® b')] ® b[d + H{b ® b% (2) 

where b®b' ~ b' ®bis the isomorphism of (classical) combinatorial R. The function 
H{b ® b') is called the energy function. We will give explicit forms of the combina- 
torial R and energy function in the next section. 

2.2 Combinatorial R and energy function 

We give explicit description of the combinatorial i?-matrix (combinatorial R for 
short) and energy function on B^'^®B^' . To begin with we define few terminologies 
about Young tableaux. Denote rows of a Young tableaux Y by yi,y2, ■ ■ - yr from 



r eib®b' a ip^ib) > e^ib') 

\ b®eib' ii^i{b) < €,,{¥), 

r f,b®b' a ^,{b) > e,{b') 

1 b®fib' iiip,{b)<e,{b'), 
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the top to bottom. Then row word row{Y) is defined by concatenating rows as 
row(Y) = yrVr-i ■ ■ - Vi- Let x = (xi, X2, . • •) and y = {yi, y2, . . .) be two partitions. 
We define concatenation of x and y by the partition {xi + yi,X2 + y2, ■ ■ ■)■ 

Proposition 2.1 ([42]) b b' e B''^' ® S^'-"' zs mapped tob'®be B''^'' ® B''^' 
under the combinatorial R, i.e., 

b®b' ®b, (3) 

if and only if 

{b' ^ row{b)) = (b^ row(b')). (4) 

Moreover, the energy function H{b ® b') is given by the number of nodes of ib' ^ 
row{b)) outside the concatenation of partitions (s^) and (s"'). □ 

We define another normahzation of the energy function H such that for b®b' & 

H{b ® b) := min(r, r') ■ min(s, s') - H{b ® b). (5) 

For special cases of i?^'* ® 5^'*', the function H is called unwinding number and H 
is called winding number in [32]. Explicit values for the case b®ib' E B^'^ ® B^'^ are 
given by 

H{b®b') = x{b<b'), H{b®b') = x{h>b'), (6) 

where x(True) = 1 and x(False) = 0. 

In order to describe the algorithm for finding b and b' from the data {b' ^ rovj{b)), 
we introduce a terminology. Let F be a tableau, and Y' be a subset of Y such that 
Y' is also a tableau. Consider the set theoretic subtraction 6 = Y\Y' . If the number 
of nodes contained in ^ is r and if the number of nodes of 9 contained in each row 
is always or 1, then 6 is called vertical r-strip. 

Given a tableau F = (6' rovu{b)), let Y' be the upper left part of Y whose 
shape is (s^)- We assign numbers from 1 to r's' for each node contained in 6' = y \y 
by the following procedure. Let 9i be the vertical r'-strip of 9 as upper as possible. 
For each node in 9i, we assign numbers 1 through r' from the bottom to top. Next 
we consider 9\9i, and find the vertical r' strip 6*2 by the same way. Continue this 
procedure until all nodes of 9 are assigned numbers up to r's'. Then we apply inverse 
bumping procedure according to the labeling of nodes in 9. Denote by Ui the integer 
which is ejected when we apply inverse bumping procedure starting from the node 
with label 1. Denote by Yi the tableau such that (Yi «— ui) = y. Next we apply 
inverse bumping procedure starting from the node of Yi labeled by 2, and obtain 
the integer U2 and tableau Y2. We do this procedure until we obtain Ur's' and Yr's'- 
Finally, we have 

b' = {(Is ^ Ur's'Ur's'-l ■ ■ ■ Ul), b = Yr's'- (7) 
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Example 2.2 Consider the following tensor product: 



b(S)b' 



e B 



2,3 



B 



3,2 



From b, we have row{b) = 235114, hence we have 





2 


3 




1 


3 


4 


^ 235114 




4 


5 





1 


1 


3 


43 


2 


2 


5 




36 


32 




45 


4i 




54 







Here subscripts of each node indicate the order of inverse bumping procedure. For 
example, we start from the node 4i and obtain 



1 


2 


3 


4 


\ 


1 


1 


3 


4 




1 


2 


3 


43 


2 


3 


5 




2 


2 


5 




2 


3 


5 




3 


4 




^ 1 = 


= 3 


3 




, therefore, Yl = 


36 


42 




4 








4 


4 






45 






5 






/ 


5 








54 







Ui 



Next we start from the node 42 of Yi. Continuing in this way, we obtain ^6^5 ■ ■ - Ui 

-. Since (0 ^ 321421) = 



321421 and Y, 
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3 


3 


4 


4 


5 


5 



we obtain 



1 


1 


4 


2 


3 


5 



2 


3 




1 


1 


3 


4 




2 


2 


4 


5 




3 


4 



3 


3 


4 


4 


5 


5 



1 


1 


4 


2 


3 


5 



2 


3 


3 


4 


4 


5 



3. 



Note that the energy function is derived from the concatenation of shapes of b and 

□ 



b', i.e.. 



3 Combinatorics on the set of paths 
3.1 Combinatorics 

3.1.1 Transportation matrices and tabloids 

Let rj, be a positive integer, a = (ai, ■ ■ ■ , a„) and /5 = ■ ■ ■ ./?„) be two compo- 
sitions of the same size. Denote by M„(a,/?) the set of matrices M = ('^i,j)i<j,j<n 
such that 

(8) 



m 



Remind that a tabloid of shape a and weight /3 is a filling of the shape a by the 
numbers 1, 2, ■ ■ ■ , n in such a way that the number of i's appearing in the filling is 
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equal to It is clear that the number of tabloids of shape a and weight (3 is equal 
to the multinomial coefficient 

A row (column) weakly strict tabloid of shape a and weight /3 is a filling of the 
shape a by numbers 1, 2, ■ ■ ■ ,n such that the numbers along each row (column) are 
weakly increasing and /3i is equal to the number of z's appearing in the filling. 

Example 3.1 Take a = (2, 1, 3, 1), /3 = (1, 3, 0, 2, 1), then 



2 4 



is a row weakly strict tabloid of shape a and weight /3. 



□ 



(m 



ijj G Mn{a, (3), we fill the row of 
. For example, let 



We denote by Tab (a, /?) the set of all row weakly strict tabloids of shape a and 
weight (3. It is well-known that there exists a bijection between the sets M„,(q;,/5) 
and Tab{a,l3). Namely, given a matrix m 
the shape a by the numbers l*"*!, 2"^'^, ■ ■ ■ 

/ 
1 





1 l\ 


2 10 
10 



Then the corresponding row weakly strict tabloid is 



4 


5 




1 






2 


2 


4 


2 





. To each tabloid T, 



one can associate the reading word, namely the word obtained by reading the filling 
of tabloid T from the right to the left starting from the top row. For example, for 
the tabloid T displayed above, w{T) = 5414222. 

If weight yu of a tabloid T is a partition, we define the charge c{T) of tabloid T 
to be the charge c{w{T)) of the reading word. See page 242 of [31] for the definition 
of the Lascoux-Schiitzenberger charge 

Example 3.2 Take standard tabloid 



T 



4 


5 




1 






3 


6 


7 


2 





then 



and therefore c{T) = 3 + 2 



wiT) = 5342lo74633i2o 
+ 4 + 3 + 1 + = 13. 



□ 
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3.1.2 Plane partitions 

Let A be a partition. A plane partition of shape A is a tabloid tt of sliape A such 
that the numbers in each row and each column are weakly decreasing. For example, 



7 


5 


4 


7 


4 


4 


3 


3 




3 







is a plane partitions of shape (3, 3, 2, 1). 

A plane partition vr has a three-dimensional diagram, consisting of the points 
{i,j,k) with integer coordinates such that {i,j) G A and 1 < A; < 7r(i,j), where 
i^{iij) G TT is the number that is located in the box (z,i) G A. By definition, the size 
|7r| of a plane partition vr of shape A is |7r| = ^ihj)- 

Let a and (3 be two compositions of the same size. Denote by VV{a,l3) the set 
of plane partitions vr such that 

^ n{i,i + k) = ^aj, ^ n{i + k,i) = ^ Pj. 

k>0 j>k k>0 j>k 

Finally, let us remind two classical results 

(A) (P. MacMahon, see e.g. [21], page 81) Let l,m,n be three positive integers, 
and B be the box with side-lengths l.m, n. Then 

y,H= TT l-ct-'-''-\ 

^ ii 1 _ Qi+j+k-2 

nCB (i,j,k)eB ^ 

(B) (Robinson-Schensted-Knuth, see e.g. [35], Chapter 3) There are bijections 

M{a, (3) VV{a, (3) Tab(a, (3). 



3.2 Paths 

3.2.1 B^'^ type paths 

Let a be a partition of size n. A path p of type B^'^ and weight a = (ai, a2, - ■ ■ , On) 
is a sequence of positive integers aia2 ■ ■ ■ an such that ctj = = i}- We denote 

by V{a) the set of all paths of type B^'^ and weight a. A path p is called a 
highest weight path if the sequence aia2 ■ ■ - On satisfies the Yamanouchi condition. 
We denote by the set of all B^'^ type highest paths with weight a. It is well 

known that the total number of B^'^ type paths of weight a = (ai, 02, ■ ■ ■ , ctn) is 
equal to the multinomial coefficient 

{ai + a2 + h On)! 

ai\a2\---aj. 



8 



and there are bijections 



Via] 



1:1 
< > 



Tab(a,r). 



Let us describe the general prescription to get the corresponding tabloid from a 
given path. Let the path aia2---a„ G 'P(a), we recursively add letters to the 
tabloid according to follows. Starting from the empty tabloid, 

assume that we have done up to aj_i and have gotten a tabloid T^*~^-*. Then we 
add the letter i to the right of the Oj-th row of T*^*~^) and get T^'^\ For example, 
the path p = 4221343 can be related to the following transportation matrix and row 
strict tabloid 



M 



/oooiooo\ 

1 1 

1 1 

1 1 



T 



4 




2 


3 


5 


7 


1 


6 



3.2.2 General rectangular paths 



More generally, we define path to be an arbitrary element of tensor product of 
crystals (g) S''^''*^ O ■ ■ ■ O 3''^''^. Recall that for type A^n^ case, S^'^ is, as the 

set, consisting of semi-standard tableaux over alphabet {1, 2, ■ ■ - rT, + 1}, and tensor 
product of crystals B ® B' set, cartesian product of two sets B and B' . 

Crystal graph structure on the set B ® B' is is given according to [19]. Weight 



A = (Ai,--- 
given by 



A„_|_i) of a path b = bi^b2 



bL^B' 



B 



r2,S2 



\i = total number of letters i contained in tableaux B'^^'^^ , ■ 
For example, 

|1|2| r^Twi r^m 

G S^'^ ® 5^'^ ® B' 



B 



rL,SL 



(9) 



?2,2 



is a path of rectangular shape R = {{2^), (2^), (2^)), and its weight is A = (4, 3, 3, 2, 2). 
Note that the number of standard (i.e., weight of (1^)) rectangular shape R = 
{{rl'')a=i,2,- ,l} paths is equal to the generalized multinomial coefficient 



N 

Rir ' ' ^ Rl 



(10) 



where N = Yla''^a-^a, and for any diagram A, Hx{q) denotes the hook polynomial 
(see definition on page 45 of [31]) corresponding to diagram A. 



Comments. Summarizing, one has the following (equivalent) combinatorial de- 
scriptions of the set of (crystal) paths of type 5^'*^ ® B^'^^ ® • • ■ ® S^'*" and weight 
a = («!, . . . , Q!l) as the set of 
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(a) transportation matrices A4l{o(, s), 

(6) row weakly increasing tabloids T(a, s), 

(c) plane partitions VV{a, s). 
For the given path hi®b2® ■ ■ ■ ®bn E B^'^^ ® B^''^'^ ® ■ ■ ■ ® 5^'*", the corresponding 
element in T{a, s) is determined as follows. Staring from the empty tabloid, we 
recursively add letters to the tabloid according to 61, 62, ■ ■ ■ , as follows. Assume 
we have done up to and have gotten the tabloid T^''~^\ Denote the number of 
k contained in bi by Xk. Then, for all k, we add letters i for Xk times to the right of 
the k-th row of T(*-i) and get T^'\ 

Example 3.3 Consider the path 



3 


3 


® 


1 


1 


2 


3 




1 


3 




2 


3 


3 



of type S^'^ ® B^''^ ® -B^'^ ® B^'^ and weight a = (3, 2, 6). The corresponding tabloid 
and transportation matrix are 





2 


2 


3 








T = 


2 


4 












1 


1 


2 


3 


4 


4 



M 




To find the plane partition which corresponds to the tabloid T (or matrix M), one 
can apply the Robinson-Schensted-Knuth algorithm [27j to the multi-permutation 



w :- 



1 1 122333333 
22324112344 



which corresponds to the tabloid T. One has 



RSK 



2 3 4 4 



2 3 3 3 



Finally, the plane partition we are looking for, can be obtained from the pair of semi- 
standard Young tableaux displayed above by gluing the Gelfand-Tsetlin patterns 
that correspond to the Young tableaux in question: 



/73IO 7310\ 
5 2 1 7 3 1 

4 2 ' 4 1 
2 3 



7. 7. 4. 3 

5. 'a 'a 'i. 

4./'2./'l..''l. 
2. '2. '1. 'g. 



3 
■5 

11 

2 '6 '8 'll 



The result is a plane partition from the set VV{{2, 4, 2, 3), (3, 2, 6, 0)). 



□ 
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Remark 3.4 For the reader's convenience, let us recall the way to get the corre- 
sponding Gelfand-Tsetlin pattern from a given semi-standard tableaux. By looking 
contents of a semi-standard tableau T, one can define a sequence of partitions 

= c /i(^) C ■ ■ ■ C = /i (11) 

such that each skew diagram fi^^' n) is a horizontal strip, see e.g.. 

Chapter I of [31]. Starting from the sequence of partitions (fTTl) . one can define the 
corresponding Gelfand-Tsetlin pattern x := x{T) by the following rule 

x^'\T) = shape{fi^'^), {l<i<n). (12) 

It is known that thus obtained x indeed satisfies the defining properties of the 
Gelfand-Tsetlin patterns. □ 

Remark 3.5 One of the basic properties of the BBS is that the second Young 
tableaiil] (of weight /3) obtained by means of the Robinson-Schensted-Knuth algo- 
rithm, is conserved under the dynamics of the BBS [1] (see also [ITJ [1] for the other 
connections between the simplest BBS and the RSK algorithm). Nowadays, con- 
served quantities and linearization parameters (or angle variables) of the BBS are 
completely determined in the most general settings [28], and, surprisingly enough, 
they are elegantly described by the so-called (unrestricted) rigged configurations 
[20} [211 [231 [261 [39l [2] ■ The latter result is a consequence of a deep theorem stated 
in Lemma 8.5 of [26]. □ 



3.3 Statistics on the set of paths 
3.3.1 Energy statistics 

For a path 6i (g) f)2 ® ■ ■ ■ ® &l e 5'"''^' ® -B''^'*'' ® ■ ■ ■ ® B''^'^^, let us define elements 
6^*^ G B^^'^^ for i < j hy the following isomorphisms of the combinatorial R; 

6i ® 62 ® ■ ■ ■ ® h-i ®bi®---® bj^i ®bj®--- 
~ 61 ® 62 ® ■ ■ ■ ® ® 6i (g) ■ ■ ■ (g) b'j'~^^ (g (g ■ ■ • 

. . . 

~ 6i(g62(g---®6i_i®6j.^^(g---(g6j._2®^j-_i® ■■■ , (13) 

where we have written bk (g b^^'^^^ ~ b^^^ (g b'/^ assuming that b^j'^ = bj. 

For a given path p = 5i ® 62 ® ■ ■ ■ ® &l e 5'''''*' ® B''^'^^ (g ■ ■ ■ ® B''^'""^, define 
statistics E{p) by 

Eip) = J2Hik^b^;-''^). (14) 

Define the statistics maj(|?) by 

maiip) =J2Hik®bf+^^). (15) 

i<j 

^ Equivalently the upper part of the corresponding plane partition. 



11 



For example, consider a path a = ai ® 02 ® ■ ■ ■ ® ol £ {B^'^)'^^. In this case, we 
have ttj*^ = Oj, since the combinatorial R act on B^'^ ® 5^'^ as identity. Therefore, 
we have 

L-l 

maj(a) = ^(L - i)xiai < Oi+i). (16) 

Define another statistics tau as follows. 

Definition 3.6 For the path p E B''^''^ ® B''^'"^ ® ■ ■ ■ ® B'''^'"^ , define r^-" hy 

r'''^(p) = maj(M(") ®p), (17) 

where uf'' is the highest element of B'^'^ . □ 

We use abbreviation r for the statistics r^'^ on B^'^ type paths a G {B^'^)^^, i.e., 

r(a) = maj(l ® a) = maj(a) + L (1 — (5i^aJ, (18) 

where ai denotes the first letter of the path a. This r is a special case of the tau 
functions for the box-ball systems [29l [36] which originate as ultradiscrete limit of 
the tau functions for the KP hierarchy |16j . 

Ej=i/^i with 
(19) 

(20) 
□ 

Note that we have a = a[i] ® ap] ® ■ ■ ■ ® a^n]? i-e., the path a is partitioned according 
to jj. It is convenient to identify as statistics on a tabloid of shape fi whose 
reading word coincides with the partitioned path according to /i. For example, to 
path p = abcdefgh and the composition /i = (3, 2, 3) one associates the tabloid 



c 


b 


a 


e 


d 




h 


9 


f 



3.3.2 Statistics charge 

Let p be a path of type B^'^, denote by T{p) the corresponding row strict tabloid and 
by w{T{p)) its reading word. Define the charge of a path p to be the charge of tabloid 
T{p), i.e. the charge of the reading word w{T{p)). For example, take p = 4221343. 
Then w{T{p)) = 4327561, and therefore, c{p) = 3 + 2 + 1+ 4 + 3 + 3 + = 16. 
If /X is a composition, define c^{p) = J2i c(p[i]), where 

cf. Definition 3.3. 
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Definition 3.7 For composition = (/zi,/i2,--- ,/in), write = 
convention /U[o] = 0. Then we define a generalization of t by 

n 
i=l 

where 



Lemma 3.8 One has 



mi + 1 



where denotes the first letter of the path p[i^ . 

Proof follows from two simple observations that 

L 



r(p) + c(l®p)= ( 2 ), c{l p) ~ c{p) = L 5i,p„ 



where L denotes the length of path p. 



□ 



Comments. It follows from [29] that on the set of semi-standard Young tableaux, 
i.e., on the set of highest weight paths, the statistics tau coincides with statistics 
cocharge. Therefore, one can consider the statistics tau as a natural extension of the 
statistics cocharge from the set of semi-standard tableaux to the set of tabloids, or 
on the set of transportation matrices. □ 



3.3.3 Haglund's statistics 

Tableaux language description For a given path a = ai ® 02 ® ■ ■ ■ ® G 

(5^'^)®^, associate tabloid t of shape fi whose reading word coincides with a. This 
correspondence is the same as those used in Definition 13. 7[ Denote the cell at the 
i-th row, j-th column (we denote the coordinate by of the tabloid t by tij. 

Attacking region of the cell at is all cells {i, k) with k < j or {i + l,k) with 
k > j. In the following diagram, gray zonal regions are the attacking regions of the 
cell (ij). 




Follow [8], define |lnvjj| by 

|lnvjj| = ^{{k,l) G attacking region for \ tki > tij}. (21) 
Then we define 

|Inv^(a)| = J2 l^^^^jl- (22) 
If we have i(t-i)j < tij, then the cell is called by descent. Then define 

Des^(a)= Yl if'^-j)- (23) 

all descent 

Note that (/Uj — j) is the arm length of the cell 
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Path language description Consider two paths a'^^\a^'^^ G [B^'^)'^^. We denote 
by a*^^^ ® a*-^'' = ai ® 02 ® ■ ■ ■ ® a-^^. Then we define 



11 1 

fc=i i=fc+i 



(24) 



For more general cases a'-^-* G (-B^'^)' 
define 

lnv(^,,^,)(a(^),a(2)) := lnv(,„^,)(a(i), 1^^"^"^^^ ® a^'^)- 
Then the above definition of |lnv^(a)| is equivalent to 



and a^^^ G (i?i'i)®/^2 satisfying /ii > //2, we 

(25) 



n-l 



|lnv^(a)| = ^lnv(^^,^^^^). 

For example, consider the following tabloid (a = 2312133212); 

|Inv(6,4)(a)| = 10. 



(26) 



3 2 



We associate the paths a^^^ = 231213 and a^^) = 1®23212. Then 

|lnv(„(i)^^(2))| =x(«i < ^2) + x(ai < 03) + x(ai < 04) + x(ai < 05) + x(«i < ^6) + 
X(«2 < 03) + x{a2 < 04) + x(a2 < as) + x(o2 < ae) + x(«2 < «?) + 
X(a3 < 04) + x(a3 < as) + x(a3 < ae) + x(a3 < «?) + x(«3 < ^8) + 
X(a4 < as) + X(a4 < Oe) + x(a4 < a?) + x(a4 < Og) + x(a4 < ^9) + 
X(as < ae) + x(«s < a?) + x(as < ^s) + x(os < Og) + x(as < aio) + 
X(a6 < ar) + x(«6 < ^s) + xi^e < Og) + x(a6 < aio) + x(«6 < an) 
= (1 + + 1 + + 1) + (0 + + + + 0) + (1 + + 1 + + 0) + 
(0 + 1 + + + 1) + (1 + + + 1 + 1) + (0 + + + + 0) = 10. 

Consider two paths a^^) G (B^'i)®^! and a^^) g (5i'^)®''2 satisfying /xi > /is. 
Denote a = a*-^-* ® a^^-* . Then define 



Ml 



Des(^i,^2)(a) = ^ {k-{fii- /is) - l)x(afc < a^+^J. 

fc=/xi-/i2 + l 

For the tableau T of shape /i corresponding to the path a, we define 

n 

Des^(T) = ^Des(^^,^^^,)(a[i] ® a^^+ij). 
Definition 3.9 ([7]) For a ^jat/i a, statistics maj^ zs defined by 

Ml 

maj^(a) = ^ maj(ti_i (g) ® ■ ■ ■ ® ^m;,*)- 



(27) 



(28) 



(29) 



i=l 
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and inv^(a) is defined by 



inv^(a) = |Inv^(a)| - Des^(a) 



(30) 
□ 



If we associate to a given path p G V{X) with the shape i^i tabloid T, we sometimes 
write maj^(p) = maj(T) and inv^(jo) = inv(r). 

Example 3.10 For highest weight paths of weight A = (2,2,2) and shape = 
(4,2), the following is a list of the corresponding tabloids associated with data in 
the form (maj^, inv^) : 



2 


2 


1 


1 


3 


3 






1 


3 


2 


1 


3 


2 





(2,4) 



:i,4) 



2 1 



(2,3) 



1 1 



(0,6) 



2 1 



:i,5) 



Let us observe that 



P6P+(A) 



which is different from 



K,,M^t) = q' + qH + qH + qH' + qH' 

Another interesting choice is A = (2, 2, 2) and /i = (3, 3). The following is a list 
of all such paths with corresponding statistics: 



2 


1 


1 


(3,3) 


2 


1 


1 


(3,2) 


1 


2 


1 


(3,2) 


3 


3 


2 


3 


2 


3 


3 


3 


2 


1 


2 


1 


(2,3) 


3 


2 


1 


(0,6) 




3 


2 


3 


3 


2 


1 



Thus we have 

^inv^(p)^maj^(p) = ^6 ^ ^3^3 ^ ^3^2 ^ 2g2^3 

peP+(A) 

which is again different from 

KxAq^ t) = (f + (ft" + 9'^' + 9'^' + g'^'- 



□ 
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4 Box-ball system 



In this section, we summarize basic facts about the box-ball system which will be 
used in the next section. For our purpose, it is convenient to express the isomorphism 
of the combinatorial R 

a®hc:ih' ®a (31) 

by the following vertex diagram: 



a- 



h' 

Successive applications of the combinatorial R is depicted by concatenating these 
vertices. 



Following [TOl [5], we define time evolution of the box-ball system T^'^\ Let 
uf* G i?"'' be the highest element and hi G B^^'^\ Here the highest element 



u 



(a) 
1,0 
,{«) 



ul ' G -B"' is the tableau whose z-th row is occupied by integers i. For example, 



u 



(3) 



1 


1 


1 


1 


2 


2 


2 


2 


3 


3 


3 


3 



Define u\'^j and b'^ G B*"''*' by the following diagram. 



bi 



u 



(a) 
LO' 



■U 



(a) 
1,1 



■u 



(a) 
1,2 ' 



U 



(a) 

LL-1' 



b' 



■ U 



(a) 
LL 



h'r 



(32) 



(a) 

by 



are usually called carrier and we set uf^ := ui"' . Then we define operator 



,(«) 



t/") ih) = b' = b[®b'^®---®b'L. (33) 

Recently EH] , operators t/""* have used to derive crystal theoretical meaning of 
the rigged configuration bijection. 

It is known ([28] Theorem 2.7) that there exists some / G Z>o such that 



T, 



(a) 



T, 



(a) 



T, 



(a) 



1+2 



(34) 



If the corresponding path is 6 G (i?^'^)®^, we have the following combinatorial 
description of the box-ball system [l5l Hi]. We regard [T] G B^'^ as an empty box of 
capacity 1, and |T] G B^'^ as a ball of label (or internal degree of freedom) i contained 
in the box. Then we have: 

Proposition 4.1 (|T0]) For a path b G {B^'^)«'^ of type aII\ T^\b) is given by 
the following procedure. 



1. Move every ball only once. 
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2. Move the leftmost ball with label n + 1 to the nearest right empty box. 



3. Move the leftmost ball with label n+1 among the rest to its nearest right empty 
box. 

4. Repeat this procedure until all of the balls with label n + 1 are moved. 

5. Do the same procedure 2-4 for the balls with label n. 

6. Repeat this procedure successively until all of the balls with label 2 are moved. 

□ 

There are extensions [TU [15] of this box and ball algorithm corresponding to gen- 
eralizations of the box-ball systems with respect to each affine Lie algebra [13l [12] . 
Using this box and ball interpretation, our statistics r(6) admits the following in- 
terpretation. 

Theorem 4.2 ([29j Theorem 7.4) For a path b G (B^^^)®^ of type A'n\ T{b) 
coincides with number of all balls 2, ■ ■ ■ , n + 1 contained in paths b, T!^\b), ■ ■ ■ , 



Example 4.3 Consider the path p = 
compute T(jj){j)) in two ways, 
(i) First we compute by Eq.( [T8|) . 

T(a) = maj(l ®a) = 7- l + 6- 
r(6) =maj(l®6) = 7 - 1 + 6-0- 



b where a = 4312111, b = 4321111. We 



-5-0 + 4- 1 
5-0 + 4-0- 



-3-0 
3-0- 



-2-0 
2-0- 



-1-0 
1 ■ = 



^11, 
7. 



Thus we obtain T(7j)(p) = r(a) + t(6) = 11 + 7 = 18. 

(ii) Next we use Theorem 14. 2[ According to Proposition 14.11 the time evolutions of 
the paths a and b follows: 



4 


3 


1 


2 


1 


1 


1 


1 


1 


4 


1 


3 


2 


1 


1 


1 


1 


4 


1 


1 


3 


1 


1 


1 


1 


4 


1 


1 


1 


1 


1 


1 


1 


4 


1 


1 


1 


1 


1 


1 


1 


4 


1 


1 


1 


1 


1 


1 


1 



4 


3 


2 


1 


1 


1 


1 


1 


1 


1 


4 


3 


2 


1 


1 


1 


1 


1 


1 


1 


4 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



Here the left and right tables correspond to a and b, respectively. Rows of left 
(resp. right) table represent a, T^\a), ■ ■ ■ , {T!^^)^{a) (resp., those for b) from top 
to bottom. Note that we omit all frames of tableaux of B^'^ and symbols for tensor 
product. Counting letters 2, 3 and 4 in each table, we have r(a) = 11, r(6) = 7 
and again we get T(^7j){p) = 11 + 7 = 18. Meanings of the above two dynamics 
corresponding to paths a and b are summarized as follows: 
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(a) Dynamics of the path a. In the first row, there are two sohtons (length two 
sohton 43 and length one soliton 2), and in the second row, there are also two 
solitons (length one soliton 4 and length two soliton 32). This is scattering 
of two solitons. After the scattering, soliton 4 propagates at velocity one and 
soliton 32 propagates at velocity two without scattering. 

(6) Dynamics of the path h. This shows free propagation of one soliton of length 
three 432 at velocity three. 

□ 

5 Main results 

5.1 Haglund— Haiman— Loehr formula 

Let H^{x]q,t) be the (integral form ) modified Macdonald polynomials where x 
stands for infinitely many variables xi,X2, - ■ ■ ■ Here H^{x\ q, t) is obtained by simple 
plethystic substitution (see, e.g., section 2 of [9]) from the original definition of the 
Macdonald polynomials [HI]. Schur function expansion of H^{x; q,t) is given by 

Hf,{x;q,t) = ^kx,i,iq,t)sxix), (35) 

X 

where Kx f^{q,t) stands for the following transformation of the Kostka-Macdonald 
polynomials: 

KxAq,t)=t<^^KxAQ,t-')- (36) 

Here we have used notation ?7,(yu) = ~ Then the celebrated Haglund- 

Haiman-Loehr (HHL) formula is as follows. 

Theorem 5.1 ([8]) Let a : ^ Z>o he the filling of the Young diagram n by 
positive integers Z>0; (ind define x^ = n^eAt '^°'(") ' -^^^'^ Macdonald polynomial 
Hfj_{x;q,t) have the following explicit formula: 

H^{x;q,t)= g'-^'^'^r^j^'^x". (37) 

o-:/j— >Z>o 

□ 

From the HHL formula, we can show the following formula. 
Proposition 5.2 For any partition yU and composition a of the same size, one has 

J2 ^-,(p)^maj,(p) ^ J2 K,,^k,,,{q, t), (38) 

where rj runs over all partitions of size \^\. 



18 



Proof. Indeed, [8], if two fillings a and a' belong to the same Soo orbit, then 
Inv((T) = Inv(a'), Des(cT) = Des(a'). ■ 

Corollary 5.3 The (modified) Macdonald polynomial H^{x\q^t) have the following 
expansion in terms of the monomial symmetric functions m\{x): 

H,{x-q,t) = Y.i Yl g'^^-^^^r^j-^^) I m,{x), (39) 
Ahl^^l \pev{x) J 

where A runs over all partitions of size |/x|. □ 

To find combinatorial interpretation of the Kostka-Macdonald polynomials K\^^{q, t) 
remains significant open problem. Among many important partial results about 
this problem, we would like to mention the following theorem also due to Haglund- 
Haiman-Loehr: 

Theorem 5.4 ([8J Proposition 9.2) // /ii < 2, we have 

kxM^t)= J2 (40) 

□ 

It is interesting to compare this formula with the formula obtained by S. Fishel [3], 
see also [21], [25] . 

Concerning validity of the formula Eq. (l40!) . we state the following conjecture. 
Conjecture 5.5 Explicit formula for the Kostka-Macdonald polynomials 

Kx,f,iq,t)= J2 g'^^^^f^r^j^W. (41) 

is valid if and only if at least one of the following two conditions is satisfied, 
(i) /Ui < 3 and /X2 < 2. 

(ii) X is a hook shape. □ 

Example 5.6 As an example, the following is the list of the tabloids associated 
with the highest weight paths of weight A = (3, 2, 1) and shape /i = (4, 2). Here we 
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also include the values of the statistics in the form (maj^,inv^) : 



2 


1 


1 


1 lO 


2 


1 


1 


1 


(1,4) 


1 


2 


1 


1 n 


3 


2 




2 


3 




3 


2 




1 


2 


1 


1 ('9 9^ 


2 


2 


1 


1 


(1,4) 


2 


2 


1 


1 n c;^ 


2 


3 




3 


1 




1 


3 




3 


2 


1 




3 


2 


1 


1 


(0,5) 


1 


1 


2 


1 (n n\ 


2 


1 




1 


2 




3 


2 




1 


1 


2 


1 l'9 n 


2 


1 


2 


1 


(1,3) 


2 


1 


2 


1 n /i^ 


2 


3 




3 


1 




1 


3 




3 


1 


2 




3 


1 


2 


1 


(1,4) 


1 


3 


2 


1 n 


2 


1 




1 


2 




2 


1 




1 


3 


2 








1 


2 





Then the generating function is 

^inv^(p)^maj^(p) _ ^6 ^ ^5^ ^ + + + qH^ + 2qH + 2qh'' + gt^ 

peP+(A) 
which is different from 

^A,M(g, t) = q^ + qH + 2g5 + ?,qH + + g^^t^ + 3^^^ + 2qH'^ + g^t + qt^ , 

Even if we consider the special value t = 1, these two polynomials are distinct. Yet 
other examples which show that the formula fHUl) does not hold if the condition (i) 
and (ii) of Conjecture 15.51 break down, is given in Example 13. 101 Let us remark that 
the choice A = (2, 2, 2) and = (3, 3) will give an example of both specializations 
g = 1 and t = 1 give distinct polynomials. □ 

5.2 Generating function of tau functions 

Our main result is an elementary proof for special case t = 1 of the formula Eq. fl38l) 
in the following form. 

Theorem 5.7 Let a he a composition and be a partition of the same size. Then, 

J2 g-^J» = 5^K,,„K,,,(g,l). (42) 

□ 

Conjecture 5.8 Let a be a composition and fi be a partition of the same size. Then, 
^-E,>.a. ^ g--(^') = ^J^,,,K,,^(g,l). (43) 

p£V{a) ri\-\fj,\ 

□ 



20 



Here X]i>2 '^i equal to the number of letters other than 1 contained in each path 
p e V{a). 

Let us remark that in view of general definition of r'"'*, our is related with 
r^'^, whereas maj^ is related with r^'^ where r is bigger than the length of weight 
a (see Section 6). As for intermediate r^''*, see Conjecture 15.191 for some further 
information. 

Example 5.9 Let us consider case a = (4, 1, 1) and = (4, 2). The following is a 
list of paths p and the corresponding value of tau function r(4 2) (p)- For example, the 
top left corner 1 111123 3| means p = |T](g)|T|®|T|(g)|T|(g)[2]®[3] and r(4,2) (p) = 3. 



111123 3 


111132 2 


111213 2 


111231 3 


111312 2 


111321 3 


112113 3 


112131 4 


112311 3 


113112 3 


113121 4 


113211 2 


121113 4 


121131 5 


121311 4 


123111 5 


131112 4 


131121 5 


131211 4 


132111 3 


211113 5 


211131 6 


211311 5 


213111 6 


231111 7 


311112 5 


311121 6 


311211 5 


312111 6 


321111 4 



Summing up, LHS of Eq. fl43l) is 

peP({4,i,i)) 



+ Aq'^ + 7q^ + 7g^ + 7g + 4. 



In the RHS of Eq.( l43|) . nontrivial contributions come from the following 4 terms: 

-^(6),{4,l,l)-^(6),(4,2)(q', ^) + -^{5,1),{4,1,1)-^{5,1),{4,2)(?,^) 

+ -f^(4,2),(4,l,l)-^(4,2),(4,2)(9i + -^(4,1,1),(4,1,1)-^(4,1,1),(4,2) (O', t) 

= l-{l) + 2-{q^ + q^ + qt + q + t) 

+ 1 ■ {2q^ + qH + q^ + 2qh + q^ + qt + t^) 

+ 1 ■ {q^ + qH + q^ + 2q^t + q^ + 2qH + qf + qt) 
= + (t + 3)g" + (3t + 4)g3 + (4t + 3)5^ + {f + At + 2) q + ^ + 2t + 1. 

By setting t = 1, we get Eq.dlS]). □ 



Proof of Theorem 15.71 

Definition 5.10 Let jj, he a composition and T he a tahloid of size \^\. Denote hy 
T^^^ the part ofT which is filled hy numhers from the interval . Then 



Lemma 5.11 



.(T) = ^c(tW). 



i>l 



Cm(T) 



(44) 
□ 

(45) 



where the sum in the right hand side runs over standard tahloids T of shape a. 
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Proof. Recall the following three formulas from [31], Chapter VI. 
Formula 1. 

K^M^t) = ir,,^,(ri,g-i)r(^)g"('^'), (46) 

and thus 

J2 K,,xK,,,{q, t) = Y, K,,xK„,it-\ q-')t<^\<^'\ (47) 
n V 
As a corollary of the formulas above, 

J2K,,xKr,AQ, 1) = ^K,,Ai^,,;.'(l,g-')g"(^'^ (48) 

Formula 2. 

J^{x-l,t) = {t,t)^>e^>{x). (49) 

Formula 3. 

it,t)Mx) = J2 ^y^.^' S,{x,t) =J2Kx,iinit)S,{x,t). (50) 

Ahr ^ Xhr 

Therefore, 

J^(x; l,t) = n E ^A(0,(iK)W'5A«(x,t) , (51) 
and after the plethystic change of variables X \ — > we obtain 



J^(x;l,t) = n E ^,w,(iMj)W^A«(x) I (52) 



Claim 5.12 

E ^.«,a.^)W^A«(x) = E^'^^^^''^'^^ (53) 

x(i)\-fir T 

where the second sum runs over all standard tabloids T of the size r, and c{T) 
denotes either the charge ofT, or the value of tau function on the path corresponding 
to tabloid. 

In the case c{T) = c{T) the charge of tabloid T, this result is due to Lascoux- 
Schiitzenberger; in the case of c{T) = n{fi') — r(T) this statement is a corollary of 
Theorem 7.4 and Corollary 6.13 from [29], where identification of tau function and 
cocharge is given. 

Corollary 5.13 

J^(x;l,t) = 5^5^tvmm,(x), (54) 
where the second sum runs over the set of standard tabloids of shape X, and c^'(T) = 

Ec(r«). 
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On the other hand, 

V 



mx{x), (55) 



and therefore 



5^X,,^(l,t)K,,A = ^tvm, (56) 

V T 

where the sum in the right hand side runs over the set of standard A-tabloids. 
Finally, 

J2 K„aK,Aq, 1) = E '^"^"'^ = E ^"^'^'^-'^^^^ (57) 

where the the third sum runs over the set of all standard a-tabloids. 
It remains to observe that according to Lemma 3.4, 



5.3 Comments on generalizations of Section 15.21 

In order to clarify nature of tau statistics, we consider possible generalizations of 
the results in Section [5^ 

5.3.1 Regularization map and parabolic Kostka polynomials 

The main objective of this Section is to give an interpretation of the energy statistics 
partition function for the BBS as the value of a certain parabolic Kostka polyno- 
mial, see e.g., [22l [25] . This observation allows to write a fermionic formula for the 
parabolic Kostka polynomials in question, see e.g. [22], as well as appears to be 
useful in the study of the BBS, see e.g. [29] . 

Definition 5.14 Letp be a path of type B"^^'^^ and weight X, define regularization 
p = reg(p) of the path p to be 

reg(p) = (1 ■ ■ ■ n - 1)^" . . . (1 . . . i)^'+i . . . l^^p, (58) 

where (1 ■ ■ ■ z) := 1 ® ■ ■ ■ C?) and we have omitted all symbols ®. □ 

Let T be semi-standard Young tableau (i.e., highest weight) corresponding to the 
regularized path p. 
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Lemma 5.15 ( [29] . Lemma 7.2) Assume that all Vi = 1, then 

t{p) = r{p) + const = c(T) + Const, 

where 



(59) 



L 



l<a<b<n 



Const = L{L - /ii) - I 2 ; ■ ^ " \ 2 

a=l ^ l<a<fe<n 

anC? L = ?77,a. □ 

Example 5.16 Let pi = 3332221, then 

Pi = 1212121113332221. 

We have r(pi) = 7, t{pi) = 46, and t{pi) — r(pi) = 39. On the other hand, let 
P2 = 3223123, then 

P2 = 1212121113223123, 

and t{p2) = 14, r(p2) = 53, so that r(p2) — '^{P2) = 39, as expected. Time evolution 
of p2 is 



3 


2 


2 


3 


1 


2 


3 


1 


1 


1 


2 


3 


1 


2 


1 


1 


1 


1 


2 


3 


1 


1 


1 


1 


1 


1 


2 


3 


1 


1 


1 


1 


1 


1 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



Corollary 5.17 



□ 



(60) 



where Ci is a constant and partitions A and k are defined as follows; A = (Ai, ■ ■ ■ , A„), 



and Ai = J2a>il^a, ^ = 1, 
Proof. Indeed, we have 



,n; K — {Ki, ■ ■ ■ , Kn), Hi — X]a>i+1' ^ ~ 



, n — 1. 



LHS = I^""^'^ = ^""^ E = Q'^'KAMiQ-'), (61) 
pe-P(A) f 

where summation in the third term runs over all Littlewood-Richardson tableaux 
of shape A and weight (k. A). ■ 

Conjecture 5.18 Let Ri = (s/')i.2,... be a sequence of rectangles, then 

J2Kr,,Rj<r,,RM = K A,i.,R,)i^~') , (62) 



where A denotes partition (Xla>j ^'^Y' > = {J2a>i+i ^'^Y' ■ 



□ 
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's well-known that parabolic Kostka polynomials K\^fi{(i), where R = {(s^"}a=i,2,. 
is a dominant (i.e., si > S2 > ■ ■ ■ ) sequence of rectangular shape partitions, satisfy 
the so-called duality theorem 

K,,n{q) = Ky,n'{q~')q''^''\ (63) 

where R' denotes a dominant rearrangement of the sequence of rectangular shape 
partitions {{s^'')}a=i,2,- 

n{R) = ^ min(ra,rfe)min(sa,Sb). (64) 

l<a<fe 

Questions 1. Can one define a set V{Ri, R2) of paths vr of rectangular type R2 = 
'S)a(^a') rectangular weight Ri = <S)b(f^b'')y energy function e(7r) such that 

KKr,,rM := J^g^W = J2k,,RiKv,rM- (65) 

□ 

Questions 2. Are there some "physical interpretations" of the duality theorems 

i^A,if(g) = g"(''^i^A',if'(^) (66) 

and 

KKn.^nM = q^^''^^ KKn^^n^iq^') (67) 

for polynomials Kx,R{q) and KKR^^R^{q)l □ 

5.3.2 Generating functions of generalized tau statistics 

Let us consider generating functions of the generalized tau statistics: 

r"'*(j9) :=maj(w(")®p), (68) 

where ui'' G B"^'^ is the highest element and p G B"^^'^^ ® ■ ■ ■ ® 5^^"^^. Let 77 be 
a partition and i? be a sequence of rectangles, denote by K^^ji{q) the correspond- 
ing parabolic Kostka polynomial. For any partition A, let Kr^^x := Kr^^xil) be the 
corresponding Kostka number. Then we have: 

Conjecture 5.19 Let R = {{ri, Si)f^i} be a sequence of rectangles. Denote by Vr 
the set of all paths (including non-highest elements) corresponding to the tensor 
product of crystals B""^^'^ ® B""^'"^ (g) ■ ■ ■ ® B''^''^. Then, 

J2 = J2 K,,R{q-')K,^x. (69) 

pePfl, wt(p)=A T? 

where s > maxj{sj}, 0^ ^ ^- D 
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In other words, the statistics r^'* defines the essentially unique class of polynomials 
although its definition depends on choice of u^s^ . As we see in the following example, 
this is not obvious from definition of r'''*. 



Example 5.20 Let us consider the case A = (4,3,3,2,2) and S^'^ ® 5^'^ ® 5^'^ 
Then we have total of 759 paths, and by direct computations, we have the following 
summation over all 759 paths. 









r+ 


8q 


^4 


-33gS 
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89g^ + 


161g^ + 


m 
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163g^ 


+ 


82g3 


+ 24g2 
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4-33g 
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f 89gi° 


+ 161g^ 


+ 19^ 
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S2g^ + 24g^ 
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f 33g 
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f 89g^ 


+ 161g^- 


F 198 




f 163g^ - 
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^ + 24g^ 
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8q 


^4 


-33gS 
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89^^ + 


161g^ + 


19J 


iq^ 
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163g^ 
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82g3 


-4- 24r/^ 


Y^qr^^^ip) 










iq'^ 


+ 


33g6- 


f J 


I9q^ + 161q^ + 198 




+ 163g2 




^2g 


+ 


24, 


where s = 2 


1 




, 5 and. 


in 


the last 


expression, 5 < 


r 


< 


10. 


However, if 


we look at 


specific paths. 


for example. 
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h 
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Then we have 
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^2,5 


-7-3,5 


-,-4,5 ^5,5 


-,-6,5 ^7,5 


-,-8,5 


-,-9,5 


-,-10,5 
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5 




7 
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3 
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2 




2 








b 


3 


7 
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2 2 
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2 
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5 




3 


3 


3 3 


3 


3 




3 







In particular, dependences of T^''^{b) on r are different for each b. □ 



Remark 5.21 We use the same notations of Conjecture 15.191 Then, as we will see 
in Corollary 16.41 below, we have 0,, = for all r > /(A), where /(A) is length of 
weight A. □ 

5.4 Generating functions related with the energy statistics 
on the set of rectangular paths 

Let us consider generating function with respect to E statistics. This is a more 
traditional problem compared with it for r^'**. Indeed, special cases of this problem 
as well as its restriction to the set of highest weight paths were considered by several 
authors, see, e.g., [32], [TTl [33], HQl HI] and references therein. 
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Conjecture 5.22 Let R = {{ri, Si)^^i} is a sequence of rectangles. Denote by Vr 
the set of all paths corresponding to tensor product of crystals B^^'"^^ ® B^^''^^ ® ■ ■ ■ ® 
B''^^'^. Then, 

J2 /^^^ = E ^^vA<l)I<v^- (70) 

pe'Pfl, wt(p)=A v^\M 

□ 

Comments. An algebra-geometric definition of tlie parabolic Kostka polynomials 
(also known as generalized Kostka polynomials) has been introduced in [13] (see 
also [25]) as a natural generalization of the well-known formula, [31], p. 244, (1), for 
the Kostka-Foulkes polynomials in terms of a g-analogue of the Kostant partition 
function. Based on the study of combinatorial properties of the algebraic Bethe 
ansatz, a fermionic formula for the parabolic Kostka polynomials has been discovered 
by the first author in the middle of 80's of the last century and has been proved in 
the full generality in [26]. A "path realization" of the Kostka-Foulkes polynomials 
has been obtained in [32], and finally, the formula 



PG^+,H>wt(p) 



has been proved in 



□ 



Example 5.23 Let us consider the case A = (4, 6, 3, 1) and i? = {(2, 2), (2, 2), (3, 2)}, 
i.e., B^'^ CS> S^'^ <E) B^'^. Then we have the following nine paths: 
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1 1 
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3 



1 


2 


2 
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Therefore, the LHS of Eq.([70D is 



Eip) 



11 



11 



12 
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On the other hand, non-zero contributions for RHS of Eq. fITOl) comes from 



V 






(6,4,3,1) 


+ 


1 


(6,5,2,1) 


^10 + ^11 


Q 


(6,4,4) 




Q 


(6,5,3) 




q + 


(6,6,2) 




+ q^ 



Summing up, 

RHS = (g^ + 5i°)-l + (?i° + g")-l + (gi°)-l + (gi^)-2 + gi2.2 
= + 3gi° + 3g" + 2g^^ 

which coincides with the LHS. □ 



6 Discussion: r'^^ and E 



So far in this paper, we have considered several statistics including generalized tau 
statistics r*"'* and more traditional one E. Let us investigate several further aspects 
of these two statistics. Our main results in this section are (i) to show that E belong 
to the class of statistics r'''* and (ii) to show that r'''* stabilize when we increase the 
value of r. The following proposition will be a key property. 

Proposition 6.1 Let bi j be the integer at the i-th row, j-th column of the tableau 
representation ofbG B"^''^' , and let the highest element of B'''^ be nJ''^ . Then we have 

H{u''' (g) 6) = (71) 

if r > br'^s' o-nd s > s' . 

Note that br'^s' is the largest integer in the tableau representation of b. 



Proof. According to the algorithm presented in Proposition l2.1l we have to compute 
the insertion b row{u^'^). This is worked out in Lemma [6.21 below, and shape of 
the resulting tableau coincides with the concatenation of two tableaux b and w*"'*. 
Hence i7(w^'* ® 6) = due to Proposition 12.11 ■ 

Lemma 6.2 Under the same assumptions of Proposition \67T\ the insertion 

b ^ rr ■■ -r ■■ -22 ■■ -211 (72) 
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gives the concatenation of tableaux u'''* and b, i.e., 



1 


1 




1 


^1,1 


&1,2 






2 


2 




2 


^2,1 


&2,2 




&2,s' 


















r' 


r' 




r' 


&r',l 


&r',2 
























r 


r 




r 



(73) 



Proof. This insertion procedure can be divided into two steps. First, we show 



b := {b ^ rr ■ ■ ■ r ■ ■ ■ 66 ■ ■ ■ 6 S5 ■ ■ ■ 5) 



bi,i 


bi,2 




bi,s' 
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6 
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b2,i 


b2,2 




b2,s' 
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6 




6 


































br',1 


br',2 




br',s' 


r 


r 


r 


r 



(74) 



where 6 = r — r' + 1 and 6 = 6+1. Note that we have 6 > 1 since (i) by the 
semi-standard property of b, we have br',s' > r' and (ii) by the assumption r > br',s', 
thus r — r' > 0. Again from the assumption r > br',s', we have 



i > br'-i x' 



(0 < Vz < r') 



(75) 



by the semi-standard property < ^fc,s'- Consider the insertion [b ^ rr- ■ - r). 

From Eg. (1751) . we have r > r — {r' — 1) > br'-(^r'-i),s' = ^i,s'- Thus the first row of 
(6 ^ rr ■■■ r) is 6i_i6i_2 ■ ■ ■ ^i.s'Tr ■ ■ ■ r, and the remaining rows are identical to the 
corresponding rows of b. If r = 1, this finishes the proof (i.e., by r > br',s' we have 
br'^s' = 1); therefore let us consider the case r > 1. Assume that we have, for some 
r > k > 6 and k = k + 1, 



b° := {b ^ rr ■ ■ ■ r ■ ■ ■ kk ■ ■ ■ k kk ■ ■ ■ k) 



bi,i 


bi,2 




bi,s' 
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b2,i 


b2,2 




b2,s' 
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k 
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b*,i 


b*,2 




b^,s 
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r 
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br',1 


br',2 




br',s' 



(76) 

where * = r' — k + 1. We insert (A; — 1) into this b°. From Eg. (1751) and the assumption 
k > 6, we have k > 6 = r — {r' — 1) > 6r'-(r'-i),s' = ^i,s', thus k — l> bi^s'- Therefore 
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inserted [k — 1) bumps k at [s' + l)-th column of b°. As the next step, we have 
to insert k to the second row of b°. By the similar reasoning, it bumps {k + 1) 
at {s' + l)-th column. In this way, insertion of {k — 1) causes downward shift of 
{s' + l)-th column of b° and addition of {k — 1) to the first row of {s' + l)-th column. 
Similarly, we see that the second insertion of {k — 1) causes shift and addition of 
{k — 1) to the (s' + 2)-th column. Continuing in this way, we see that successive 
insertions {b° ^ {k — l){k — 1) ■ ■ ■ {k — 1)) gives the tableau b° with replacement k 
hj k — 1. By induction, we can show Eg. ( 17^ . 

Next, we consider the insertion of {6 — 1)* = {6 — 1){6 — 1) ■ ■ ■ {6 — 1) into b: 

l:={b^{6-ir). (77) 

We denote the i-th row of b (resp. b) by bi (resp. bi). Although we independently 
repeat insertions of {S — I) for s times, we can argue more systematically as follows. 
Let us consider the first row of b: 

bi = bi,ib^,2---bi,s'^l^. (78) 

s 

As we saw in the last paragraph, we have bi^s' < <^ from Eq. flTSj) . Suppose there are 
ki letters 6 in the first row of b, i.e., bi^s'-ki+i = &i,s'-A;i+2 = ■ ■ ■ = bi^s' = Thus 
there are (s + ki) letters 6 in the first row of b. After inserting {6 — 1)^, precisely s 
letters S are bumped and go down to the second row, and the first row becomes 

k = &i,i&i,2 ■ ■ ■ f'l,.'-^, {6-l){6-l)---{6-ll 66---6^ . (79) 

s ki 

In particular, ki letters 6 on the right of the first row of b are precisely reproduced 
in the right part of bi. 

Now we have to consider the insertion of s letters 6 which are bumped from bi 
into the second row of b: 

h = b2,ib2,2 ■ ■ ■ b2,s' (<5 + 1)(5 + 1) ■ ■ ■ (5 + 1) . (80) 

^ v ' 

S 

Again, from Eg. ( 1751) . we have 62,5' ^ S + 1, and suppose that there are ^2 letters 
{6 + 1) in the second row of b, i.e., 62,s'-fc2+i = ^2,s'-fc2+2 = ■ ■ ■ = 62,5' = 6+1. After 
inserting 6^, precisely s letters {6 + 1) are bumped and go down to the third row, 
and the second row becomes 

k = 62,1^2,2 ■ ■ ■ b2,s'-k, 66--- H6 + l){6+l)---i6 + ll . (81) 

Again, k2 letters {5 + 1) on the right of the second row of b are precisely reproduced 
in the right part of 62- 
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As we see in the above discussions, this procedure can be continued recursively. 
Let the number of 6 + i contained in the {i + l)-th row of b be /cj+i- Then fej+i is 

bi+i,ibi+i,2 ■ ■ ■ i6 + i-l){S + t-l)---{S + t-l){5 + i){5 + t)---{6 + i) 

^ ^ ^ ' 

s ki+i 

(82) 

and the right /Cj+i letters {5 + i) are copy of those originally contained in the right 
of {i + l)-th row of b. In this way, b contains copy of the letters in b. 

As we have investigated insertion of {6 — 1)*, let us consider the insertion of 
{6 — 2Y = {6 — 2){6 — 2) ■ ■ ■ {6 — 2) into b. Consider the row bi. Suppose that there 
are li letters {6 — 1) within the first row of b, i.e., bs'^(ki+h)+i = ■ ■ ■ = b^'-fci = S — 1. 
Thus there are total of /i + s letters 5 — 1 in the row bi. Therefore, after insertion 
of (5 — 2)*, the row bi becomes 

&i,i&i,2 ■ ■ ■ bi,s'-ik,+i,) {S -2){5-2)---{5- 21(5 - 1){5 - 1) ■ ■ ■ [5 - 1)^ M ■■ ■ 5^. 

s li ki 

(83) 

In particular, li letters {6—1) and ki letters 6 are copy of the corresponding letters 
of the first row of b. As the result, s letters {6 — 1) are bumped from bi. Therefore, 
for the row 62, we have to insert (5 — 1)*, and again obtain copy of the letters 6 and 
{6 + 1) contained in the second row of b. We can continue this procedure until the 
bottom row of b. Therefore each row of the resulting tableau {b ^ (5 — 2)*) contains 
copy of at most two species of letters in the original b. We can recursively continue 
insertions of {6 — 3)*, {6 — 4)*, ■ ■ ■ , 1*, and each insertion generates copy of part of 
letters of b. Finally we get the result Eq. fl73l) . ■ 

The following theorem shows that the statistics E essentially belong to the class 
of statistics r'''*. 

Theorem 6.3 Let p = bi®b2® ■ ■ ■ ^bi e B''^'""^ (g) B''^'"^ ® ■ ■ ■ ® B''^'''^, and denote 
by N the largest integer contained in tableau representation of p. If r > N and 
s > maxjjsj}, we have 

E{p) = C-T^'%p), 

C = min (rj, rj) ■ min (sj, Sj) 

i<j 

Proof. Recall the definition of r^'* 

r^'*(p) =maj («'■'* ®j9), (86) 

where maj is defined by Eq.( IT5l) and m^'* is the highest element of S^'*. Within the 
definition of r**'*, let us first consider the energy functions involving u^'^ and next 
consider the remaining ones. As for the terms involving n'''*, we have 

H{u'^'' ®bf) = (l<Vi<L) (87) 
31 



(84) 
(85) 



by Proposition 16.11 On the other hand, for the remaining contributions, recall that 
bi G i?**^'*' and b^j^^^ G B^^'^K Then from the definition of normalizations, we have 

H{hi ® b^l^^^) = min(ri, rj) ■ mm{si, Sj) - H{bi O bf^^^). (88) 
Combining both contributions, we obtain the sought relation. ■ 

Corollary 6A Let p = bi (g) b2 (S) ■ ■ ■ (S) bi G B''^'''^ (g) B"'^'"^ ® ■■■ ® B""^^"^ , and 
denote by N the largest integer contained in tableau representation of p, and define 
S = maxj{si}. Then we have 

t''''{p)=t^'^{p) (89) 
for all r > N and s > S . □ 

Therefore Conjecture l5.19l means that we can define at most independent statistics 
r''''*(p), where N is the largest integer contained in the path p, and these statistics 
have essentially unique generating function. 

Let us remark physical interpretation of r^'''. For the paths (including non- 
highest elements) of shape 5^'"*^ ® B^''^^ ® ■ ■ ■ ® B^''^^ and statistics t^'^ {S = 
maxj{sj}), there is a straightforward generalization of Theorem 14. 2[ see Section 
4.1 of [29]. Under the same assumptions, t^'^ is identified with cocharge of the un- 
restricted rigged configurations. It will be an interesting problem to find a physical 
interpretation of the more general T^''^{p). 
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